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ABSTRACT 

I study the relationship between the spatial distribution of intergalactic metals and 
the masses and ejection energies of the sources that produced them. Over a wide 
range of models, metal enrichment is dominated by the smallest efficient sources, as the 
enriched volume scales roughly as E 3 ^ 5 ~ M 3 / 5 while the number density of sources 
goes as ~ M . In all cases, the earliest sources have the biggest impact, because fixed 
comoving distances correspond to smaller physical distances at higher redshifts. This 
means that most of the enriched volume is found around rare peaks, and intergalactic 
metals are naturally highly clustered. Furthermore, this clustering is so strong as to lead 
to a large overlap between individual bubbles. Thus the typical radius of enriched z ~ 3 
regions should be interpreted as a constraint on groupings of sources rather than the 
ejection radius of a typical source. Similarly, the clustering of enriched regions should 
be taken as a measurement of source bias rather than mass. 



Subject headings: intergalactic medium - galaxies: evolution 



1. Introduction 

It is now clear that the intergalactic metals detected at z ~ 3 have only a secondary impact 
on further structure formation. Numerous studies of the C iv , Si iv , O vi and C in content of 
the intergalactic medium (IGM) have been carried out (e.g. Songaila & Cowie 1996; Aracil et al. 
2004; Aguirre et al. 2005) only to find that this material is a small fraction of the metals produced 
(Pettini 1999); there have been many detailed measurements of z ^ 5 galaxy outflows (e.g. Pettini 
2001; Schwartz &: Martin 2004), but the intergalactic metal distribution is observed to be roughly 
constant over this entire range (Songaila 2001; Pettini et al. 2003); and metal ejection has been 
incorporated into several simulations (e.g. Thacker et al. 2002; Springel & Hernquist 2003; Cen et 
al. 2005), which found that it has a negligible impact on IGM cooling and the statistical properties 
of the Lyman-alpha forest (Theuns et al. 2002; Bruscoli 2003). In short, the IGM metals that we 
see are not doing very much. 
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Yet as a tracer of the higher-redshift interplay between galaxies and the IGM, intergalactic 
metals are unparalleled. It is indeed remarkable that this material has made its way from the 
centers of stars into the lowest-density environments detectable (Schaye et al. 2003; Aracil et al. 
2004), with far-reaching implications. The depletion of metals through outflows directly impacts the 
galaxy mass-metallicity relation (e.g. Dekel & Woo 2003; Tremonti et al. 2004); outflows suppress 
the formation of nearby objects (Scannapieco et al. 2000; Sigward et al. 2005); and the distribution 
of IGM metals is closely linked to the evolution of the first generation of stars (e.g. Bromm 2003; 
Scannapieco et al. 2003). 

Extracting the details of each of these processes from observation, however, requires an un- 
certain extrapolation from z ~ 3 to much higher redshifts. One approach to this problem is to 
focus on deriving constraints from the IGM composition, which can be related to the star forma- 
tion history, initial mass function, and metallicity of the sources (e.g. Aguirre et al. 2004; Qian & 
Wasserburg 2005). In this case the primary complications are due to uncertainties in abundances 
and the ionizing background. 

A second method relies on the spatial distribution of intergalactic metals. Regardless of the 
details of IGM enrichment, it is clear that metals were formed in the densest regions of space, regions 
that are far more clustered than the overall matter distribution. Furthermore, this "geometrical 
biasing" is a systematic function of mass and redshift (e.g. Kaiser 1984), and thus the observed 
large-scale clustering of metal absorbers encodes information on the scales of the objects from which 
they were ejected. Likewise, the size of each enriched region is dependent on the energy at which 
the metals were dispersed. 

Thus recent measurements of the sizes of enriched regions (e.g. Rauch et al. 2001), the metal- 
licity as a function of IGM density (e.g. Schaye et al. 2003), the C iv absorber correlation function 
(Rauch et al. 1996; Pichon et al. 2003), and the galaxy-C iv cross-correlation function (e.g. Adel- 
berger 2003), are already providing useful constraints on simulations of metal enrichment. Yet 
such detailed comparisons are too expensive to carry out over a large range of parameter space, 
and provide little intuition. One is left wondering if perhaps there might be some general rules of 
thumb that can be kept in mind when interpreting observations or selecting simulation parameters. 

It is this issue that I take up in this Letter. Adopting a simplified model, I show that the sizes 
and clustering of the enriched regions that we see should not be interpreted as ejection radii of 
individual sources or correlated with masses at z ~ 3, but that nevertheless these quantities have 
natural interpretations related to the properties of the sources of IGM enrichment. Throughout, 
I assume a cold dark matter cosmological model with parameters h = 0.65, Slo = 0.3, S7a = 0.7, 
fib = 0.05, as = 0.87, and n = 1, where h is the Hubble constant in units of 100 km/s/Mpc, Qq, 
f&A> an d &b are the total matter, vacuum, and baryonic densities in units of the critical density, <t| 
is the linear variance on the 8/i _1 Mpc scale, and n is the "tilt" of the primordial power spectrum 
(e.g. Spergel et al. 2003), with the transfer function taken from Eisenstein & Hu (1999). 

The structure of this work is as follows. In §2 I describe a general model of cosmological 
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enrichment and apply it in §3 to derive the clustering properties and sizes of enriched regions under 
a wide range of assumptions. I conclude with a short discussion in §4. 

2. Modeling Cosmological Enrichment 

As we are interested in making general statements, I adopt here an extremely simple approach. 
All outflows are taken to be pressure-driven spherical shells, expanding into the Hubble flow (e.g. 
Ostriker & McKee 1988). To keep this model as transparent as possible, I do not attempt to 
include secondary effects, such as cooling, a stochastic star formation rate, external pressure, or 
the gravitational drag from the source halo (Madau et al. 2001; Scannapieco et al. 2002). In this 
case the evolutionary equations are 

R s = —^---^-{R s — HR s ) 2 — VL m — — — , 
ptt s Ms * 

E b = L(t) - 4irR 2 s R s P b , (1) 

where the overdots represent time derivatives, the subscripts s and b indicate shell and bubble 
quantities respectively, R s is the physical radius of the shell, E b is the internal energy of the hot 
bubble gas, P b is the pressure of this gas, and p is the mean IGM background density. I assume 
adiabatic expansion with an index 7 = 5/3 such that Pi, = E b /2irRl. When outflows slow down 
to the sound speed, the shell is likely to fragment. At this point I let the region expand with the 
Hubble flow. 

In this model, the evolution of the shell is completely determined by the mechanical luminosity, 

L(t) = imL & f,f w NM b Q(t m -t), (2) 

where /* is the fraction of gas converted into stars, f w is the fraction of the supernova (SN) kinetic 
energy that is channeled into the galaxy outflow, AT is the number of SNe per solar mass of stars 
formed (each assumed to explode with 10 51 ergs of kinetic energy), M b is the baryonic mass of 
the galaxy in units of solar mass, and isN = 5 x 10 7 years. Following Scannapieco et al. (2003), I 
define the product, /*/ u ,£kin.A/', as the "energy input per unit gas mass" £ g . Assuming /* = 0.1, 
f w = 0.1, and 1 SN per 300 solar masses gives a fiducial estimate for £ g of 10 -4 ' 5 , although I vary 
this parameter over a wide range below. Note that the choice of isN has no direct impact on the 
results. 

By combining eqs. (1) and (2) with the standard analytical mass distribution, one can compute 
the porosity, which is defined as the product of the number density of sources and the bubble volume 
around each source: 

where V(£ g , M, z, z') = Airr s (£ g , M', z, z') 3 /3 and r s (£ g} M', z, z') is the comoving radius at a red- 
shift 2 of a shell from a source with total mass M' that hosts an outflow at a redshift z' with 
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an energy input per unit gas mass £ g . Finally, d ^I^ z , is the differential Press-Schechter comoving 
number density of objects forming as a function of mass and redshift calculated from: 

dn p ,,2 /9 diner 2 

= ve ' (4) 

(MM (2vr) 1 /2M cflnM V ; 

where v(M, z) = 6 c /[D(z)a(M)], D(z) is the linear growth factor, a(M) 2 is the variance associated 
with the mass-scale M, and 5 C = 1.69. Note that although, strictly speaking, ^wrrp accounts for 
both the creation of new sources and the destruction of older sources by merging into larger objects 
{e.g. Benson et al. 2005), it is sufficiently close to the formation rate for the objects in which we 
are most interested here. 

The porosity Q can be thought of as the average number of outflows impacting a random 
point in space, and it depends on only two free parameters: £ g and the minimum mass, M m \ n {z). 
Here I assume for the fiducial case that efficient star formation occurs only in halos with virial 
temperatures above 10 4 iT, because smaller objects are photevaporated after reionization, and before 
reionization, gas cooling in these objects requires H2, which is an inefficient coolant (Madau et al. 
2001) that is easily dissociated (Haiman et al. 1997; Ciardi et al. 2000). In our cosmology, this gives 
M m - m = 2.4 x 10 7 [(1 + z)/10]~ 3//2 M o , although I also consider variations in M min in §3.2 below. 

Carrying out similar integrals as in eq. (3) we can compute "porosity-weighted" estimates of 

the properties of typical sources. For example the source mass that contributes most significantly 

to enrichment can be estimated as 

f°° f°° dn 

(log(M)) Q = Q' 1 J m dM' J z dz 1 ^-^y(£ g ,MU/)log(M'). (5) 

Similar averages can be used to compute the typical comoving bubble radius, {h)q, the number 
density of sources, ^^3?) > an ^ tne source °i as (&)q> where b(z,M) = 1+ [i/(z',M r ) 2 — l] /S c (Mo 
& White 1996). 



3. Results 

3.1. Variations in Energy Input 

The properties of three enrichment models with T v i r > 10 4 K, and widely varying £ g values 
are shown in Figure 1. In all cases, the average source mass closely follows the minimum allowed 
value. The reason for this behavior can be seen from a simple Sedov- Taylor estimate. In this case 
the physical radius goes as R cx (E/p) 1 / 5 , where E is the energy of the blast and p is the ambient 
density, and thus the comoving volume goes as 

y 6 ocr 6 3 ocM 3 / 5 (l + z) 6 / 5 . (6) 

From eq. (4) the comoving number density goes as 1/M (with only a small correction from the 
term) such that Q cx M~ 2 / 5 . Notice that this is a general result, which follows from dimensional 
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Fig. 1. — Properties of T V - 1T > 10 4 K enrichment for an energetic model with E g 
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(left 



column), the fiducial model with f g = 10 (center), and a weak outflow model with £ g = 
10 -5.5 (right). Top row: Porosity-weighted average source mass as a function of redshift (solid) 
as compared to the minimum mass (dotted). Second row: Porosity- weighed average bias (b)n 
(dashed), and bias normalized versus z = 3, x D(z)/D(3) (solid) which tracks the evolution 
of the correlation function as in eq. (7). Third Row: Volume- weighed average comoving outflow 
size, (r s )q (dashed), and the typical comoving size of enriched region, max |r over i ap , (r s )g| (solid), 
which can be significantly larger due to overlapping sources. Bottom Row: Simple 1 — exp(— Q) 
estimate of the filling factor (solid), and 1/2 of the collapse fraction (dotted), an estimate of the 
enriched gas falling back onto further generations of sources. 
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analysis. The smallest sources naturally dominate the enrichment process, an effect that is only 
amplified by complications such as the additional gravitational drag in large halos (Scannapieco et 
al. 2002) or the inefficiency of OB associations at driving winds from large galaxies (Ferrara et al. 
2002). In fact only an extremely strong increase in efficiency with mass (<? g oc M n , with n > 2/3) 
can alter this conclusion. 

Similarly, the strong redshift scaling in eq. (6) means that for any given mass the earliest sources 
enrich most efficiently. This is due to the simple fact that a fixed comoving distance corresponds 
to a smaller physical distance at higher redshift, yet it has profound implications for the resulting 
spatial distribution. In the second row of Fig. 1, I plot the porosity-averaged bias for each of the £ g 
models considered. The strong increase of Vf, at high-redshift means that the relatively rare, high-z^ 
sources contribute most strongly to Q. This results in a high bias, as such rare sources are highly 
clustered relative to the lower-f peaks that collapse at lower redshifts (e.g. Mo & White 1996). 
Again this is a general result, that arises from dimensional arguments, and the strong clustering 
seen in this figure can only be amplified by IGM transitions such as reionization, which remove 
lower-mass, lower-redshift (and hence lower-i/) peaks (e.g. Klypin et al. 1999). Note that in general 
while bias increases with redshift, the amplitude of the correlation function is given by 



(b) Q D(z)]\ (r), (7) 



where £o( r ) is the matter correlation function, linearly extrapolated to the present, and D(z) 
decreases with redshift. Thus £(r, z, (b)n) remains roughly constant with redshift. 

The strong clustering of enriched regions must be also taken into consideration when inter- 
preting the typical sizes of enriched regions. In the third row of Figure 1, I plot (r s )g, which again 

1/5 

scales roughly as £ g , and is always smaller than 300 comoving kpc. In order to estimate the 
sizes of typical measured regions, however, we must also consider the overlapping of such bubbles. 
This can be estimated by considering the distance from the center of a typical source at which 
the product of the number density of neighboring sources n(r over i ap ) and the volume around each 
source is equal to one. In our formalism this gives 

1 + f (^overlap, Z, (&l)q)] X Q = 1, (8) 

where here it is more appropriate to calculate bias in the Lagrangian coordinate system that does 
not include the peculiar motions of the sources, as these motions were not included in our outflow 
model. This means (b L ) Q = (b) Q - 1 (e.g. Mo & White 1996). 

Despite the large range in £ g values considered, the overlap radius is well above {t s )q at z < 6 
for all three models. Thus the observed radii of enriched regions correspond to the sizes of typical 
overlapping groupings of sources, rather than the ejection radii from individual objects. Similarly, 
the growth of this scale with time is not caused by the expansion of material from a typical source, 
but rather is due to the formation of ever-larg er groupings of overlapping bubbles. Indeed ?"ovcriap 
increases drastically at late times in all three models, greatly outpacing the growth of the typical 
outflow as defined by {t s )q . 
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redshift 

Fig. 2. — Properties of enrichment models as a function of minimum mass. Each row is labeled 
by the assumed M m { n value, while the curves in each row are as in Figure 1. In the left column 
M m m = 10 9 Mq, and I consider two cases with £ g = 10 -4 ' 5 and £ g = 10~ 3 ' 5 . In the center panel, 
■Mmin = 1O 11 M0, and I take even higher energies of £ g = 1CP 3 ' 5 and £ g = 1CP 2 ' 5 , which still give 
very low filling factors. Finally the right column shows the results of an extreme low-mass model 
in which T vir > 10 3 K, and with low £ g values of 1CT 4 - 5 and 1CT 5 ' 5 . 
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The last row of Fig. 1 provides simple estimates of the filling factor of enriched regions. If 
the distribution of sources were completely uncorrelated, this could be computed directly from the 
porosity as F = 1 — exp(— Q), which corresponds to the solid lines in these panels. However, the 
strong overlapping between bubbles means that the true filling factor is probably somewhat smaller 
than this estimate. Finally, the dotted lines in these panels are 1/2 of the total collapse fraction, 
which is intended as a rough estimate of the gas that would fall back onto new sources in a more 
detailed simulation. This shows that infall is not important in these models, except for perhaps a 
small correction in the £ g = 10 -55 case. 

3.2. Variations in Minimum Mass 

The second major parameter in our models is the minimum mass, which is raised to 10 9 M G 
in the cases shown in the left column of Figure 2, as might be caused by the increase in the IGM 
temperature following reionization, for example. This naturally raises (log(M))g by over an order 
of magnitude, but the bias is much less affected, shifting from ~ 3.5 up to ~ 4.0. Thus although 
enrichment now occurs later, it is dominated by sources with similar v values as in the T vir > 10 4 K 
case. Also as in the fiducial models, r over i ap > r s in the observed redshift range, such that the scale 
of enriched regions is set by groupings of objects. Interestingly, as the filling factors are smaller in 
this case for the same values of £ g , raising the minimum mass actually lowers the sizes of typical 
enriched regions. However, there is probably a significant infall correction in this model. 

These effects are intensified in the extreme 1O 11 M0 model shown in the center column of Figure 
2. In this case, the bias rises to ~ 5.5, but the filling factors are so low that r ovcr iap < r s . The 
observed bubble radii are even smaller than in the 10 9 M Q case, but now one is actually looking 
around individual sources. However, the very high energies and low filling factors make this case 
unlikely. Infall only makes this worse. 

Finally, the right column of this figure shows the results of an extreme low-mass model in 
which the minimum virial temperature has been reduced to 10 3 K. Even in this case (6}q is about 
3.5 at z = 3. The radii of enriched regions are again set by clustering, and are even larger than in 
the fiducial T v ; r > 10 4 K cases with the same choices of £ g . 

4. Discussion 

It is clear that the simplified models described above gloss over many of the detailed issues that 
are now beginning to be addressed numerically. Nevertheless, this simplicity serves to highlight 
how many counterintuitive observational trends can be understood from general arguments, which 
can be explored in more detail with future simulations. 

Thus, the seeming contradiction between widespread outflows from large z ~ 3 galaxies and the 
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lack of evolution in C iv number densities is most likely related to the V oc E 3 / 5 scaling of outflows 
and the 1/M scaling in the number density of sources. Likewise the strong clustering of metal line 
systems is likely to be be reconciled with the efficient ejection of metals from small M <^ 10 10 ' 5 Mq 
galaxies (Tremoni et al. 2004) because cosmological enrichment is dominated by early, rare sources, 
that can expand to cover the same comoving volume in shorter times. Finally, as strong clustering 
results in a significant overlap between sources, this may account for the large observed sizes of 
typical enriched regions, which require extremely large ejection energies to be explained by single 
sources (Kollmeier et al. 2003). Although there is much more to be learned, it is clear that the 
efficiency of lower mass, high-redshift sources and the overlap between bubbles should always be 
kept in mind when interpreting measurements of the spatial distribution of intergalactic metals. 

I thank Andrea Ferrara, Pavel Kovtun, Crystal Martin, Piero Madau, Michael Rauch, and an 
anonymous referee for helpful comments and useful conversations. This work was supported by the 
National Science Foundation under grant PHY99-07949. 
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